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ABSTRACT 

We have applied the Tremaine- Weinberg method to 10 late-type barred spiral galaxies using data 
cubes, in Ha emission, from the FaNTOmM and GHaFAS Fabry-Perot spectrometers. We have combined 
the derived bar (and/or spiral) pattern speeds with angular frequency plots to measure the corotation 
radii for the bars in these galaxies. We base our results on a combination of this method with 
a morphological analysis designed to estimate the corotation radius to bar-length ratio using two 
independent techniques on archival near infrared images, and although we are aware of the limitation 
of the application of the Tremaine- Weinberg method using Ha observations, we find consistently 
excellent agreement between bar and spiral arm parameters derived using different methods. In 
general, the corotation radius, measured using the Tremaine- Weinberg method, is closely related to 
the bar length, measured independently from photometry and consistent with previous studies. Our 
corotation/bar-length ratios and pattern speed values are in good agreement with general results 
from numerical simulations of bars. In systems with identified secondary bars, we measure higher 
Ha velocity dispersion in the circumnuclear regions, whereas in all the other galaxies, we detect flat 
velocity dispersion profiles. In the galaxies where the bar is almost purely stellar, Ha measurements 
arc missing, and the Tremaine- Weinberg method yields the pattern speeds of the spiral arms. The 
excellent agreement between the Tremaine- Weinberg method results and the morphological analysis 
and bar parameters in numerical simulations, suggests that although the Ha emitting gas does not 
obey the continuity equation, it can be used to derive the bar pattern speed. In addition, we have 
analysed the Ha velocity dispersion maps to investigate signatures of secular evolution of the bars in 
these galaxies. The increased central velocity dispersion in the galaxies with secondary bars suggest 
that the formation of inner bars or disks may be a necessary step in the formation of bulges in late-type 
spiral galaxies. 

Subject headings: galaxies: spiral - galaxies: kinematics and dynamics - galaxies: ISM - galaxies: 
structure - galaxies: evolution - galaxies: individual (IC 342, NGC 2403, NGC 4294, 
NGC4519, NGC 5371, NGC 5921, NGC 5964, NGC 6946, NGC 7479, NGC 7741) 



1. INTRODUCTION 



The theory of resonan t structure in d isk galaxies in 
its original linear form by Lindblad ( 1963 1 or in its more 



* The derived Ha rotation curves for all galaxies are made 
available in electronic format at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr or via http://cdsweb.u-strasbg.fr/ 



complete form by Lin & Shu ( 1966 1 implies that a density 
wave pattern in the stellar disk acts on the rotating gas to 
form stars continually along spiral shock lines associated 
with the arms or virtually straight shocks associated with 
a bar. Additionally, resonances play an important role in 
the mechanism of maintaining density waves, which are 
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the back-bone of the spiral structure in the inner parts 



as well as outer regions of spiral galaxies (Lin|1970 1. One 



corollary is that over significant ranges of galactocentric 
radius the angular velocity of the density wave pattern 
may well be virtually invariant. These invariant patterns 
give rise to the bar (or each bar if nested bars are present) 
and to the spiral arms, and their respective angular ve- 
locities are termed pattern speeds, fl p . 

Observational derivation of the f2 p , together with the 
conventional rotation curve, allows us to characterise the 
main dyn a mical parameters of the disk. Tremaine & 



Weinberg (1984) suggested a virtually model indepen- 



dent purely kinematic way to derive fl v of which asim- 



plificd form was given by Merrifield & Kuijken (19951, 
where the Tremaine- Weinberg method can be written as 

n pS m(i) = {v)/(x), 

where i is the inclination angle of the disk and the pa- 
rameters (V) and (X) are the luminosity-weighted aver- 
ages of the projected velocity at each pixel, and the pro- 
jected distance from the minor-axis, respectively, mea- 
sured along a slit parallel to the line of nodes. The 
method is based on three assumptions, and only under 
these assumptions can the relation between (V) and (X) 
be used to determine f2 p : 

1. The disk is flat. 

2. The disk has a well-defined pattern speed. 

3. The surface brightness of the tracer obeys the con- 
tinuity equation. 

To make use of the Tremaine- Weinberg method, we 
need to take a spectrum along such a slit which gives 
us the projected velocity, and also to take an image giv- 
ing the luminos ity in the component used to measure th e 
spectrum (e.g., Merrifield fc Kuijken|1995]|Corsini]2 008). 
In practice, obtaining several slit spectra of galaxies is a 
tedious process, and consequently, two-dimensional ve- 
locity fields can be used to simulate a large number of 
slits, enhancing the signal-to-noise ratio and better ex- 
ploring the regime of multiple £l p s. To fulfill the third 
assumption, observations yielding stellar kinematics are 
considered ideal. However, high resolution observations 
of the stellar line-of-sight velocities over ent ire galax 



ies are still technically challenging (sec, e.g., Emsellem 



20071. Alternatively, emission-lines can be used to de- 



rive the kinematics of the interstellar medium in a much 
easier way. However, the interstellar gas is subject to 
large scale shocks and smaller scale expansions powered 
by OB associations, and hence does not fulfill the third 
assumption. 

Several studies have shown that the Tremaine- 
Weinberg method can be applies to CO velocity fields 
( [Rand fc W allin 2 004] |Zimmer et al.||2004|) or Ha veloc- 
ity fields ( Hernandez ct al. 2005b; Emsell em et al.||2006| 
2007bl |Chemin k Hernandez||2009|l^ab-| 
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m agreement with numerical simulations 



2009 1 in individual galaxies to derive f2 p s 

Furthermore, 



these fipS could be used to identify the location of reso- 
nances where predicted morphological features, such as 
bar-spiral transition or rings, have been identified. The 
latter studies have suggested that although Ha does not 
satisfy the third assumption, it may come close enough. 



Here we have applied the Tremaine- Weinberg method 
on Fabry-Perot data cubes in Ha emission which pro- 
vide, over the full extent of the disks, complete Ha kine- 
matic maps, red continuum maps, as well as continuum- 
subtracted Ha surface brightness maps for a sample of 10 
late-type barred spiral galaxies. Although we are aware 
of the limitations of such an analysis, the wide spatial 
coverage combined with the high spatial and spectral 
resolution makes this data set excellent for this type 
of analysis. When combined with morphological anal- 
ysis, this yields an angular speed parameter for the bar 
which we will show to be a good approach to the pattern 
speed. Support fo r this method can be found in Hernan- 



dez et al. (2005b I who found a main £} P for NGC4321 in 
goo d agreement with CO studies by Wada et al.| fl!998|) 
and|Rand fc Wa llin (20041. Furthermore, Emsellem et al.| 
( |2006| showed that thederived fi p for NGC 1068 is fully 
consistent with the stellar velocity field. 

In section [2] we describe the observations, followed by 
the analysis of the Ha kinematic maps to derive rotation 
curves and angular frequencies (section k3| . In section H 
we outline how we apply the Tremaine- Weinberg method 
and consider the position angle uncertainties and disk 
cove rage limitations deta iled by Rand & Wallin (2004) 



and Zimmer et al. 



(20041. The results are presented in 
section [6] compared with morphological analysis of in- 
frared images for our sample galaxies (section [5| , and 
with pre vious results from numerical simulations (sec- 



tion 



tion|7r 



5.3 1. Finally, the main conclusions are given in sec- 



2. OBSERVATIONS 



The observations were taken with the Fabry-Perot in- 
terferometer of New Techn ology for the Observat oire du 
mont Megantic (FaNTOmM, |Hernandez et al.|200"3| at the 
3.6m Canada-France-Hawaii Telescope (CFHT) and the 
1.6m Observatoire du mont Megantic (OmM) (see Ta- 
ble fl]). In addition, IC342 was obse rved with Galaxy 



Ha Fabry-Perot System (GHaFAS, |Hernandez et al. 
20081 mounted on the 4.2m William Herschel telescope 
on La Palma. All details regarding the FaNTOmM obser- 
vation conditions and instrument set up have been pub- 
lished in three papers by the Laboratoire d'astrophysique 
expeerimentale - LAE - at Universit y of Mon treal ( Her- 
nandez et al. | [2005a |Chemin et al. 2006| |Daigle et al 
2006ap . The GHaFAS data for IC 342 were obtained 
with a mean seeing of 1.7", on January 16th 2008, dur- 
ing 3200 sec integrating 16 cycles of 40 channels each, 
with an interference order (at Ha) p — 765 ± 1, mean 
finesse Ti — 15.1, and a free spectral range of 8.623 A. 

To generate a homogeneous set of Ha kinematic maps, 
and to gain more field coverage, we re-red uced the data 
using t he IDL data reduction package of |Daigle et al.| 
(2006b I, wi th consistenc y check s using the procedure 
outlined 
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hereafter Paper I). 



^ ( [2007a 

All data cubes were reduced and spatially re-binned to 
a minimum amplitude-over-noise ratio of 7. Enough 
frames were obtained for each galaxy to reach the back- 
ground red continuum on either side of the Ha line. 
In Fig. [I] we illustrate for all the galaxies, maps of the 
continuum-subtracted Ha intensity, Ha velocity, and Ha 
velocity dispersion (ct) maps over the entire field-of-view 
as indicated in Table [T] For IC342, the GHaFAS data 
only cover the central « 2' radius, and only the FaNTOmM 
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Fig. 1. — The complete data sets for the 10 late-type spiral galaxies with B-band digital sky survey images included for comparison. The 
continuum-subtracted He* emission-line intensity, line-of-sight velocity Vi os , and velocity dispersion <r are found in the top row, and in the 
bottom row we show the galactocentric radial profiles of velocity dispersion, the Tremaine- Weinberg method, the rotation curve derived 
by means of tilted-ring and harmonic decomposition of the observed velocity field, and the angular frequency analysis. In all maps, north 
is up, and east is to the left. In NGC4519 and NGC6946, we mark the points from the region inside the ILR in blue colour to illustrate 
that they lie on a different slope than the other points. The solid ellipse outlines the 25th magnitude radius of the galactic disk, the 
dotted line, the r(CR) of the bar or oval distortion, and the black solid ellipse oriented along the outer disk position angle, when present, 
outlines the position of the ILR assuming the outer disk projection parameters. The dashed line over-plotted on each velocity field outlines 
the kinematic major-axis <j>k and the solid line outlines the photometric major-axis <f> p retrieved from the RC3 catalogue. Note that the 
GHoFaS data for IC342 only cover the central 2' radius, and therefore we only display the FaNTOmM maps here. 
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Fig. 1. — continued.. 



maps are displayed here. Corresponding S-band images 
from the Digitized Sky Survey (DSS) are also presented 
for comparison with the Fabry-Perot maps. 

3. DERIVING ROTATION CURVES AND VELOCITY 
DISPERSION PROFILES 

We quantify the observed velocity fields by using the 
tilted-ring method comb ined with the harmonic decom- 

Schoenmakers et al.|1997 Fathi 

"TTie 
to 



position form alism (e.g. 



et al.||2005[ ). We assume that 
dominant feature and that 
positions on a single inclined disk, and divide the veloc- 



circular rotation is 
our measurements refer 



ity field in concentric rings and fit the systemic velocity 
V sys and disk geometry (position of the centre (xo,yo), 
inclination i, line of nodes position angle <j)). We fit the 
circular and non-circular motions within each ring simul- 
taneously. The inclinat ion is fixed at the value t aken 
from the RC3 catalogue ( de Vaucouleurs et al.|1991 1 and 



the projected position of the dynamical centre, V sya , and 
kine matic line of no des ((/>& ) are fitted iteratively (see 
also |Begeman 1987| ) . After each iteration, one of the 
parameters is fixed to the average value of all rings. 
We first fix the dynamical centre to the average xq 
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and yo value for all rings, second, we fit Vi os — V sys + 
sin(z) [ci cos((/)fc) + Si sin(0^)] for each ring. After calcu- 
lating (j>k the average value of the line of nodes for all 
rings, we derive the rotational and radial velocity com- 
ponents together with the higher harmonic terms simul- 
taneously (up to and including the third-order term). 
In the case when higher terms are absent, the harmonic 
parameters are c\ and s± which include V ro t and VJ-ad- 
This procedure minimises the contribution from sym- 
metric streaming motions caused b y bars or spiral arms 
on the derived rotation curves (see Schoenmakers et al. 



1997 Fathi et al.|2005 and Paper I, for a more details). 

To investigate how the choice of data parameters in 
the data reduction procedure affects the final results, we 
applied this procedure to velocity fields with minimum 
amplitude-over-noise values varying between 1 and 15. 
We find that the differences in the results from these data 
cubes could be used to derive reasonable error estimates, 
since the kinematic parameters derived from these differ- 
ent maps vary within a range larger than the statistical 
errors derived for each map (see Table [1. The disk ge- 
ometry that best fits each observed velocity field is then 
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Fig. 1. — continued. 



used to derive V^-ot and to section the a maps into the 
corresponding tilted rings. For each ring, we average all 
values from individual pixels to construct the a profiles 
presented in Fig. [I] 

Our kinematically derived V^yg values are similar to 
the systemic velocities listed in the RC3 catalogue (see 
Fig. I2J, however, since almost all galaxies are subject 
to significant amounts of non-cosmological velocities, we 
cannot use the V sys to directly derive cosmological dis- 
tances. For all galaxies, we have been able to find, in 
the literature, accurate distance measurements from ra- 



dio source studies, Cepheid variables, the Tully-Fisher 
relation, or c orrected for the Virgo clust e r perturbation 
(|TullylpJ88l iFreedman k, Madorel|1988l [Gavazzi et al. 
|1999| |I)ale et al.||2000[ |Kennicutt eHjrf|2003p . We use 
these distances throughout our analysis. 

3.1. Position Angle of the Outer Disk 

When applying the Tremaine- Weinberg method, the 
slits must be placed along the line of nodes, which can 
be approximated by the kinematic and/or photometric 
major- axes. 
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TABLE 1 

Observed galaxies in the sample with coordinates, magnitudes, and systemic velocities from the RC3 catalogue 



Object 


R.A. (J2000) 


Dec. (J2000) 


Type 


B-band 


V B ys 


Telescope 


Pixel Size 


Field of View 




h 


: m 


: s 


d : 


m : 


s 




mag 


1 —1 

km s 




"/pix 


' X ' 


IC 342 


03 : 


46 : 


49.7 


+68 


: 05 


: 45 


.SXT6.. 


9.10 


-4 


OmM/WHl 


1.6/0.2 


13.7 X 13.7/3.4 X 3.4 


NGC 2403 


07 : 


36 : 


54.5 


+65 


: 35 


: 58 


.SXS6.. 


9.12 


107 


OmM 


1.6 


13.7 x 13.7 


NGC 4294 


12 : 


21 : 


17.5 


+11 


: 30 


: 40 


.SBS6.. 


12.47 


421 


CFHT 


0.48 


4.1 x 4.1 


NGC 4519 


12 : 


33 : 


30.6 


+08 


: 39 


: 16 


.SBT7.. 


12.21 


1180 


CFHT 


0.48 


4.1 x 4.1 


NGC 5371 


13 : 


55 : 


40.7 


+40 


: 27 


: 44 


.SXT4.. 


11.34 


2575 


CFHT 


0.48 


4.1 x 4.1 


NGC 5921 


15 : 


21 : 


56.5 


+05 


: 04 


: 11 


.SBR4.. 


11.66 


1457 


CFHT 


0.48 


4.1 X 4.1 


NGC 5964 


15 : 


37 : 


36.4 


+05 


: 58 


: 28 


.SBT7.. 


12.60 


1382 


CFHT 


0.48 


4.1 X 4.1 


NGC 6946 


20 : 


34 : 


52.1 


+60 


: 09 


: 15 


.SXT6.. 


9.64 


7 


OmM 


1.6 


13.7 x 13.7 


NGC 7479 


23 : 


04 : 


57.2 


+12 


19 


: 18 


.SBS5.. 


11.47 


2394 


CFHT 


0.48 


4.1 x 4.1 


NGC 7741 


23 : 


43 : 


54.1 


+26 


: 04 


: 32 


.SBS6.. 


11.66 


755 


CFHT 


0.48 


4.1 x 4.1 



Listed are also the telescope on which FaNTOmM or CHoTaS were mounted, the pixel size, and the field of view for each object. Although 
the effective field of view of GHoFaS is 3.4 X 3.4 arcminu tes, due to rotation of t he observed field, we have used only the central pa 2' field 

jHernandez et al.|2008[ l. 
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Fig. 2. — Comparison between outer disk line-of-node position 
angle derived from photometry and Ha kinematics (top) and V ByB 
from the RC3 catalogue and those derived from our Ha kinematics. 
Error bars for the kincmatically derived parameters are illustrated 
and the diagonal line shows the 1:1 correspondence. 

The orientation of the kinematic major-axis derived 
from the quantification scheme described in section [3] is 
the result of the second iteration in the procedure for 
deriving V m t, i.e., after determining and fixing the po- 
sition of the dynamical centre (xo,yo) and V sya . Sim- 
ilar to the standard tilted-ring method, the algorithm 
minimises the expression x 2 = — ( c i + s i)' which 
can be seen ideally in cases when Vi os is the result of 
a circular cosinusoidal term (ci) and a radial sinusoidal 
term (si). In the presence of weak non-circular motions 
caused by weak perturbations or large beam smearing ef- 
fects as in radio observations, this approximation is rea- 
sonable. However, for our data, where local velocities 
from strongly star forming regions as well as significant 



TABLE 2 

Uncertainties for the dynamical centre position 
(xo (err), yo (err), EAST AND NORTH-WARD respectively), 
INCLINATION, (j>k VALUES GIVEN IN NORTH-EAST DIRECTION, AND 
DISTANCES AND ARCSEC-TO-PARSEC SCALES AS DESCRIBED IN 
SECTION [3] 



Object 


xg(err) 


2/0 (err) 


incl. 


<t>k 


Distance 


Scale 




it 


11 


deg 


deg 


Mpc 


pc/" 


IC 342 


1.3 


3.3 


25 


61 + 8 


3.9 


19 


NGC 2403 


2.8 


1.9 


60 


136 + 8 


3.2 


16 


NGC 4294 


1.5 


1.1 


65 


160 + 3 


14.0 


68 


NGC 4519 


0.5 


1.6 


36 


167 + 3 


15.1 


73 


NGC 5371 


1.8 


1.1 


48 


+ 12 


37.8 


183 


NGC 5921 


1.6 


5.7 


15 


147 + 9 


25.2 


122 


NGC 5964 


1.2 


2.1 


20 


109 + 9 


24.7 


120 


NGC 6946 


1.0 


1.3 


2d 


246 + 3 


6.5 


32 


NGC 7479 


1.1 


1.2 


42 


24+7 


32.4 


157 


NGC 7741 


2.6 


0.3 


48 


2 + 6 


12.3 


60 



flows along the bars in the objects could significantly 
skew the isovelocity contou r s towards th e bar axis (e.g. 



Duval 1977 Huntley 1978 Buta 19871, we are bound 



to derive misplacements between the kinematic line of 
nodes position angle <pk and the photometric counter- 
part 4> p . One step forward in treating this problem is the 
expansion of the minimisation routine to the expression 
X 2 = V[o S - s m=i( c m + s m)- However, patchiness of ob- 
served field, the presence of strong dust content internal 
to the observed galaxy, and asymmetric drift cause con- 
siderable computational complic ations which are bey ond 
the scope of this paper (e.g., 



Weijmans et al. 2008} 



We derive the outer disk line of nodes position angle 
(cj)p ) fr om the near infrared images presented in sec- 
tion IBTT] by using the average position angle for the outer 
isophotes over a 2-4 kpc range in radius. We compare 
the values with those obtained from the RC3 catalogue as 
well as previous studies of our sample galaxies, and find 
that in most cases, the photometrically derived values 
are in fair agreement (see Table [3l, however, for IC342 



NGC 5921, and NGC 6946, the <pp values from the liter- 
ature are based on comprehensive analysis, we use these 
angles for deriving the pattern speeds. 

Furthermore, we compare the x 2 values for the lin- 
ear fits to the (V), (A") pairs, and as it is expected that 
the Tremaine- Weinberg method delivers the least scatter 
along the correct line of nodes, we use the x 2 values to as- 
sess the (j) p we have decided to use. Fig. [3] illustrates the 
reduced \ 2 values for a range of position angles, where 



Pattern Speeds of Bars and Spirals From Ha Velocity Fields 



9 



IC342 



NGC 2403 



NGC 4294 



NGC 4519 



NGC 5371 




117 143 
<t> [deg] 



170 




117 137 
i [deg] 
NGC 5964 



157 



115 



135 155 
(j) [deg] 



175 




125 145 165 
[deg] 
NGC 6946 



185 



6 






20 








15 


4 














10 


2 






5 



211 





139 



159 179 199 

<t> [deg] 
NGC 7479 




15 35 
<t> [deg] 




160 180 200 
i [deg] 



Fig. 3. — Reduced x 2 values for the linear fits to the Tremaine- Weinberg diagrams along different line of nodes. The vertical dotted line 
indicated the angle <f> where the minimum x 2 value is obtained, and the dashed (red) line shows the best estimate for disk tf> p as described 
in section|3.l| For all V[ os measurements, a constant velocity error of 10 km s _1 has been used to calculate the reduced %2 values. 



IC342 



NGC 2403 




10 



30 50 
* [deg] 



NGC 5921 



70 




NGC 4294 




117 137 157 
i [deg] 

NGC 5964 




NGC 4519 



185 




NGC 5371 



50 
40 
30 
20 
10 





139 



NGC 6946 



159 179 
* [deg] 

NGC 7479 



199 



-22 



170 



50 
I 40 
£ 30 
| 20 

A 

v 10 b*H 
> of 

115 



a 50 
| 40 

£ 30 
I 20 

A 



135 155 
* [deg] 



175 



10 


211 




231 251 
i> [deg] 



271 




-2 18 38 
4> [deg] 

NGC 7741 



15 35 
[deg] 




140 



160 180 
4> [deg] 



200 
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show 



the minimum \ 2 values correspond to position angles al- 
most identical to or to within 10 degrees from <p p . This 
figure also shows that extracting ((V), (X)) pairs along 
offset position angles of ±10 degrees still delivers reason- 
ably linear fits to the (V), (X) diagrams such as those 
shown in Fig. [I] 

We use these uncertainties to measure differences in 
the £l p by changing the orientation of the nominal slits. 
The (V) /(X) j sin(i) values illustrated in Fig.Hshow that 
for (f> p uncertainties of ±10 degrees, we can bracket the 
fl p values to within the limits shown listed in Table |4| 
Allowing the slit widths to vary between 1 and 10 pixels 
implies additional differences in the derived £l p of the 
order of 10%. 

Previous analysis of the f2 p uncertainties due to incor- 
rect sl it ori entations have been car ried out by[Dcbattista 
( |2003[ ) and |Rand fc Wnllin| pUOT). who found that « 5 
degrees position angle (PA) offset is sufficient to cause 



J7p values that differ by up to 30%. In Fig. [5] illustrated 
as error bars, we show the difference in the measured £l p 
between an estimate with the correct line of nodes, and 
an estimate using an incorrect line of nodes, displaced 
by 10 degrees, as a function of the number of slits used 
in these estimates. The figure shows that this difference 
grows quickly as the number of slits falls. When the 
number of slits has fallen to around 15, matching this 
condition, the differences have grown to 30%, a n d thu s 
consistent with the findings of Rand & Wallin (20041 



However, we note that in our derivations of the \l p we 
use over 100 slits per galaxy. 

4. APPLYING THE Tremaine- Weinberg METHOD 

For each galaxy our two dimensional velocity field al- 
lows us to select a set of pseudo-slits which we can simply 
call "slits". We select these slits such that they are three 
pixels wide and parallel to the line of nodes of the outer 
disk. For each independent point along a slit we de- 



10 



Fat hi et al. 



TABLE 3 

Disk line of node position angles from different sources given in degrees from north to east, and the 4> p value which we 

use throughout the paper. 



Galaxy 



(RC3) 4>k 4> (NIR) 0(xLJ <h 



</>(literature) 



IC342 
NGC 2403 
NGC4294 
NGC 4519 
NGC 5371 
NGC 5921 
NGC 5964 
NGC 6946 
NGC 7479 
NGC 7741 



127 
155 
169 

8 

130 
145 

25 
170 



61 
136 
160 
167 


117 
109 
221 
24 
182 



48 
123 
153 
148 

11 

97 
111 
231 

22 
158 



42.5 
124.5 
147.5 
158.5 
8.0 
130.0 
142.5 
241.0 

17.5 
157.5 



40 
123 
153 
148 

11 
130 
141 
241 

22 
158 



39, 40, 87 (Newton 



1980 



124 ( Barbera et al.|2004| 

: 48 jPaturel et aTpOOOl iDistefano et al. 



Crosthwaitc 2002 



Buta Sz McCall 



1999) 



151, 

152, 145 (IPaturel et al.|2000l iRubin et al.|l999 ' 



4, 17 I Paturel et al 



149 ( Springob et al 



150, 130 Paturel et 



240 ( Springob et al. 



29, 22 (IPaturel et al.|2000 
162 (IPaturel et al. 12000 ~ 




Barbera et al. 



2001) 
Knapen et al.|2006 I 



1990) 



Laine et al.ll998 i 



termine the luminosity-weighted value of V\ os - V^ ys and 
the luminosity- weighted position (X), using as well the 
continuum-subtracted Ha or the stellar continuum sur- 
face brightness maps from the data cube as our weighting 
function. Checking that there are measurable velocity 
data from more than 50 pixels along any given slit, we 
then compute the Trcmainc- Weinberg integral along that 
slit, giving us a single point in our plot (see Fig.Q. The 
resulting plots have far more points in them than con- 
ventional plots obtained using the long-slit spectrum of 
the stellar component or CO observations, because we 
are able to use far more slits; a typical number for the 
galaxies analysed in this article is > 100. As we can see 
from Fig. [T] "clean" linear fits are found for almost all 
the objects. 

Even though all the Ha-emitting gas does not satisfy 
assumption 3, we note two points in support of using the 
Tremaine- Weinberg method with our data: 

Firstly, from the maps displayed in continuum- 
subtracted Ha in Fig. [I] (as well as those discussed in 
Paper I), a major fraction of the pixels (> 80%) in our 
maps are due to the diffuse ionised gas component, which 
is much less contaminated by the local dynamical effects 
such as effects of supernovae and winds from in dividual 
H H region s. Deep Ha observations of NGC 7793 (Dicaire 
et al.|2008| , and imaging of the SINGS sample ( (Jey et al. 



further confirm the contribution of the diffuse com- 
ponent. Combined with the typical star formation rates, 
< 1 M q yr -1 for half of our sample, estim ated by (Ken- 



nicutt et al . 2003), Martin fc Friedri| ( |1997j ), and |GohzaIez 
Uelgado et al.| ( |1997| , we expect that only a minor frac- 
tion of the velocities in the observed fields are affected by 
local effects such as winds from star forming complexes. 

Secondly, the gaseous motions which do not share the 
velocity field of the global pattern may, in large part, can- 
cel due to global symmetries (arms, bar) or local sym- 
metry (the expansion of individual H II regions) . The 
bi-symmetric streaming motions along the bar and the 
arms must cancel quite well in order to restrict the mean 
gas inflow rates towards the centre of a galaxy. Local in- 
terruptions to continuous flow, while real, are on average 
circularly symmetric in the plane of the galaxy and any 
departure from regularity which might cross out conti- 
nuity will cancel because there is sufficient symmetry in 
these departures. Example for such departures include 
star formation at the ends of a bar or two regular arms 



littered with Hn regions. Thus, the resulting values for 
the pattern speed are not necessarily less accurate than 
those obtained using the stellar component, especially as 
the latter are subject to relatively poor signal-to-noise 
conditions, and would need exposures more than one or- 
der of magnitude longer to give similar signal-to-noise 
ratios to those we find here. 

Regarding any possible effects of the spiral density 
wave in triggering star formation, there are two points 
to note. (I) Using the stellar continuum as our weighting 
function (see Table [4] and further discussion on this point 
in section |6| systematic effects from Ha due to the pres- 
ence of newly formed stars are ruled out. Furthermore, 
when we use the Ha emission intensity for weighting the 
velocities, columns 3 and 4 in Table [4] show that the re- 
sulting fi p do not differ too greatly from those obtained 
with the red continuum measured from the Fabry-Perot 
data cubes. The typical difference is 3 km s _1 kpc -1 
which fits within the typical uncertainties for our mea- 
surements. This suggests that star formation triggering 
effects, although locally important, do not change the 
overall weights applied to the velocities. (II) We consider 
velocity deviations due to density waves by making sure 
that we derive consistent values when we move along a 
given slit or across the disk. This strongly suggests that 
such velocity deviations do cancel to first order. 

4.1. Changing the limits and the weights when 
calculating the Tremaine- Weinberg integrals 

The ideal form of the Tremaine- Weinberg method uses 
the luminosity-weighted velocities and the distances in 
the full radial range of the galaxy disk, i.e. the integrals 
range from — oo to oo. Observations, on the other hand, 
cover only a limited radial range, and ours reach typically 
the r25 radius of the galaxy disk. It is thus important 
to see whether our observations are deep enough so that 
they do not yield incorrect fl p . 

We note that the Tremaine- Weinberg method gives 
pattern speeds strictly when the limits of the integrals 
reach well beyond the bar. However, in this test, we re- 
strict the radial range of integration in progressive steps, 
starting from an innermost radius of 10 pixels and in- 
creasing in steps of 10 pixels until the entire observed 
disk is covered (see Fig. pi. Each derivation yields a dif- 
ferent (V)J(X) value, treated as if it were a val ue for £l v 
(see Fig. m) . As shown in the detailed analysis by Zimmer 
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TABLE 4 



Pattern speeds (in km s kpc~ ) derived using continuum-subtracted Ha maps and continuum maps. The values in the 

BRACKETS LIST THE DIFFERENCES IN THE DERIVED PATTERN SPEEDS WHEN VARYING THE PSEUDQ SLIT ORIENTATIONS BY ±10 DEGREES, AND 

THEY ARE OF THE SAME ORDER FOR THE COLUMNS 2 AND 5. 



Obiect 


f2„ along (t>h Cl-o 


along 4> v (Ha) 


Q„ X sin(j) along <p v 


Cl v along 4>v (continuum) 


fin Literature 


— rn — 

(1) 


(2) 


(3) 


(4) 


Tel 

(5) 


7a\ 

(6) 


IC 342 


36 


31 (-5. 5) 


13 


32 


10-40 




22 


22 (-1, 9) 


19 


25 


10.5 


NGC 4294 


39 


43 (-12, 2) 


39 


37 




NGC 4519 


13 


20 (-8, 6) 


12 


33 




NGC 5371 


16 


14 (-1, 5) 


10 


17 




NOT 5921 


13 


13 (-2, 2) 


9 


10 


15-20 


NGC 5964 


23 


25 (-5, 1) 


9 


24 




NGC 6946 


21 


25 (-6, 6) 


10 


27 


22-66 


NGC 7479 


18 


18 (-3, 3) 


12 


15 


27-100 


NGC 7741 


33 


19 (-6, 8) 


14 


18 


30-145 


Columns 2, 3, 


and 4: il p derived along 


4>k and <pp using 


continuum-subtracted Ha 


emission-line intensities, and 


inclination-corrected 



(column 6) are further discussed in section [672] 



et al. (20041 these values decrease asymptotically as the 
disk coverage increases. In particular the differences be- 
tween the outermost values of these trials and the values 
obtained using the complete disks are much smaller than 
those between the values obtained limiting the field to 
a few tens of pixels and those found when using the full 
disk region. Furthermore, we confirm that by using all 
the pixels inside the slit (i.e., without the radial trunca- 
tion as shown in Fig. [6j|, the derived f2 p s are much more 
stable . This is in support of the result by |Zimmer et al.| 
(20041 who found that most of the variation in the de- 
rived pattern speed occurs within an inner radius which 
they set to 100" for their observations. 

A pattern speed derived over a radial range signifi- 
cantly shorter than the bar length will not give a mean- 
ingful value since the integrals do not converge, so we do 
not wish to justify here any variations (or indeed absence 
of variations ) calculated w ithin these ranges. Moreover, 
as found by |Shetty et al. (20071 for M51, quasi-steady 
state continuity is not satisfied by a single valued pat- 
tern speed, so it is important to analyse the derived f2 p 
variations to make sure they converge on a unique value, 
which can then be taken as the mai n value for a singl e 
and dominant m = 2 perturbation ( Henry et al.j 2003 1 . 
In our work we have used integrations over the full ob- 
served field for all galaxies to derive the pattern speeds 
presented in Tableland Fig. [I] 

In an ongoing study of the Tremaine- Weinberg method 
based on N-body + SPH numerical simulations carried 
out b y Isabel Perez (Perez et al. in preparation, see 
also Beckman et al. 20081 we find that in the case 
of one constant pattern speed in a disk, scans of un- 
truncated pseudo-slits in the inner regions as well as 
outer regions, all recover the input fl p of the simulations. 
When truncating the slits to cut the outer sections of the 
disk (compare with Fig. Rjl, the Tremaine- Weinberg in- 
tegrals applied to inner truncated slits result in different 
(V) /(X) values when compared with using the full simu- 
lated galaxy. This is a clear demonstration of the neces- 
sity of using full slit coverage and deep observations (even 
in the inner regions) due to the non-converging integrals. 
However, as expected, we have found that as soon as the 



radial disk section is comparable with or larger than the 
size of the bar, the £l p input into the simulations can be 
recovered. This feature is shared by real data as well as 
by numerical simulations. 

Furthermore, we have applied the Tremaine- Weinberg 
method using the simulated velocity field weighted by 
the simulated density maps as well as the simulated 
star-formation maps and found that although using the 
density maps provide more accurate results, the star- 
formation maps also deliver comparable (to within the er- 
rors) pattern speeds, implying that the Ha velocities are 
reliable when applying the Tremaine- Weinberg method. 
For the galaxies presented here, Table [3] shows that the 
Ha and continuum deliver comparable £l p (note that the 
continuum maps used for this exercise also come from 
our data cubes). The exception is NGC 4519, where the 
continuum image assigns higher weights to the velocities 
inside the bar, and thereby introduces significant differ- 
ence between the fl p derived using the continuum and 
Ha maps. The higher pattern speed of the bar, using 
the continuum image delivers a higher total f2 p . 

In the presence of multiple structures with different 
pattern speeds, when deriving the (V)/(X) for an in- 
ner component, such as an inner bar inside a large bar, 
by applying the notional slits on the inner parts of a 
two-dimensional map, the majority of the pixels will be- 
long to the disk region outside the bar. These pixels will 
therefore affect the derived pattern speed of the inner 
structure. This effect was tested by e.g., IZimmer et al. 



(2004 their Fig. 7 ) who found that once the slit cov- 
erage is comparable with the size of the main bar, the 
Tremaine- Weinberg method delivers reliable f2 p s. More- 
over, cutting the disk in radial bins seems to reveal the 
presence of multiple pattern speeds when the disk sec- 
tion is comparable to structure with a correspondingly 
distinct pattern speed. In Fig. [T] we illustrate the (V), 
(X) points from the inner regions of two galaxies where 
this effect can be seen. However, as the derivation of 
secondary pattern speeds is more complicated (see e.g., 
Meidt et al.|2008 l, this test is to be used with caution in 



those cases. 

Finally, we confirm that the £l p differences using the 
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Fig. 5. — Comparison between Q p and the varying width of the 
pseudo slits. For a comparison between our work and the CO work, 
the vertical dotted line indicates when a total of S3 15 slits covers 
the entire observed field (which is the typical number of slits used 
in the previous studies). Lowering the resolution (or using wider 
slits) delivers Op unc ertainties in agreement with the findings of 
pebattista] ( |2003[ > and |Rand fc Wallin| f2004| l . 
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Fig. 6. — Illustrating the case of IC342 as an example: Sections 
of the disk used to explore the radial dependency of the Tremaine- 
Weinberg method to investigate whether the depth of the obser- 
vations could deliver erroneous values. The diagonal lines indicate 
the orientation of the fiducial pseudo-slits (one in every ten slits 
is marked), and the cross in the bottom left panel shows the posi- 
tion of the dynamical centre. Written on each panel, are the radial 
extent and number of extracted slits. 

continuum images instead of the Ha surface brightness 
maps, does not significantly change the location of the 
resonances throughout our analysis. 

5. COMPARISON WITH MORPHOLOGY 

Morphologically, the bar radius r(bar) is commonly 
determined by analysing ellipse fittin g profiles (e.g., 
Athanassoula| 1992 Wozniak et al.|1995[ ), or the two- fold 
symmetric com ponent in a Fourier decomposition of th e 
isophotes ( Elmegreen fc Elmegreen|l985b| |Butafl986a| . 
We have used both methods to verify the validity of the 
Corotation Radius r(CR) values that we have derived 
from the Tremaine- Weinberg method applied to the Ha 
velocity fields. 

5.1. Ellipticity Profiles 



Empirical studies by Erwin ( 2005 ) a nd simulations 



by |Martinez-Valpuesta et al. 
( |20U7 1 have shown that in systems 



pdnEf and jGadotti et 
systems with single bi 



al. 
bars. 

careful examination of ellipticity (e = 1 — b/a) profiles 
can deliver a reliable m = 2 bar radius. For our sam- 
ple galaxies, we analyse near infrared if s -band im ages 
from the 2MASS catalogue (|Jarrett et al.||2000) and 
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Fig. 7. — The variation of the (V)/(X)/ sin(i) values, along the tf> p angle, using different sections of the disk as described in section |4,l| 
The diamond in the IC342 panel shows the value derived using the inner 2' field of GHoTaS, and the horizontal solid line indicates tne 
pattern speed derived using the entire observed field with upper and lower limits indicated by the dashed lines (column 3 of Table [1J. 



deep 3.6 /im images from the Spitzer archives (PLFazio, 
PLKennicutt, PI: Kenney) in order to derive the bar 
length. For the galaxies where Spitzer data were not 
available, we retrieved deep ii-band image from the 2.2m 
Calar Alto telescope (NGC4519; PI: Boselli) and deep 
if -band images from the 4.2m William Hcrsch el Tel e- 
scope (NGC5921 and NGC5964; |Knapen et aTp003| |. 

We detect the presence of a bar within a galaxy disk 
from the ellipticity and line of nodes position angle <t> pro - 
files traced by the galaxy isophotes ( Jedrzejewski||1987[ | . 
The transition from the (round) bulge-dominated centre 
to the disk is typically characterised by an increase in 
ellipticity Following Martinez- Valpuesta et al. (2006 1, 
we define the semi-major axis length of the bar (r{bar) e ) 
as the radius where the ellipticity value has fallen by 
15% (empirically derived from the nu merical simulations 



As an alternative to analysing the ellipticity profiles, 
the isophotal shapes can be de scribed by means of the 
Fourier decomposition method ([E lmegreen & Elmegreen 



by Martinez- Valpuesta et al. 2006J ) from the radially 



1985b| |Buta||1986alb||Aguerri et~aT11998p . Certain pref- 
erence can be attributed to this method as it avoids con- 
fusion caused by the morphological changes that accom- 
pany the onset of spiral structure in the disk, which could 
artificially lengthen the bar. Simulations predict that 
Fourier decomposition of the azimuthal luminosity pro- 
file is dominated by the ( m = 2) at the ba r and falls 
sharply when the bar ends ( |Ohta et al.| [l990), however, 
in reality, matters are more complicated, and a combina- 
tion of the Fourier components or the phase shift of the 
seco nd component can be used to derive the length of a 
bar dAthanassoula fc Misiriotis|2002||Aguerri et al.|2000 ' 
[20031 j Zhang k Buta)|2007| ~ 



outermost local maximum accompanied by no significant 
change in <f> (see Fig.|8|. It is important that no position 
angle variation is detected around this radius, as twists 
could imply erroneous ellipticity diagnostics. 

We compare the r(CR) derived from kinematic analy- 
sis with the bar radius derived from the images (see Ta- 
ble [5]), and find that the ratio r(CR)/r(bar) e a grees with 



We apply the Fourier decomposition to the deprojected 
Spitzer, William Hcrschel, and Calar Alto near infrared 
images described in section 5.1[ and 
and forth components (I 2 , 



in Fig. 



resent the second 
8] (see 



Elmegreen 



the predictions f rom numerical simulation s ( Athanas 
soula et al.||1990| lO'Neill k, Dubinski||2003[ ). This result 



further suggests that peak ellipticities could underesti- 
mate the radial extent of the bar, though not by a large 
factor, since the ellipticity profiles shown in Fig. [8] illus- 
trate that the radius where ellipticity value has fallen by 
15% is usually not located far from the radius of peak el- 
lipticity. Furthermore, for almost all galaxies, we derived 
comparable ellipticity profiles using the 2MASS -fG-band 
images. However, these images are not deep enough to 
deliver good quality ellipticity profiles at large radii. 



& Elmegreen 1985b, for a detailed presentation of the 
method). In this paper, since we are mainly interested in 
the lengths of the bars in our sample, we do not analyse 
the Fourier decomposition extensively, but only use it to 
ensure that these profiles indicate the presence of a bar 
in agreement with those derived from the ellipse fitting 
analysis. We use a combination of the following criteria 
to determine the bar lengths (r{bar)p), and where these 
three criteria agree, we can specify the end of the bar 
quite well. 

1. The bar ends after the m = 2 component has been 
dominant, and starts decreasing accompanied by 
an increase of the odd components. This indicates 
that the two-fold symmetry of the bar component 
has decreased. 



5.2. Fourier Decomposition 



2. The c ontrast parameter defined by |Ohta et al.| 
fll990], hi la = (l + h + h + h) / (l-h + h-h) 
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TABLE 5 

COROTATION RADIUS r(CR) AND BAR LENGTHS DERIVED FORM ELLIPSE FITTING (r(6ar) e ) AND FOURIER ANALYSIS OF THE NEAR INFRARED 
IMAGES WITH UNCERTAINTIES CALCULATED BY ALLOWING FOR THE fl p TO VARY WITHIN THE LIMITS LISTED IN TaBLe|3] In NGC4519, 
r(CR) FALLS OUTSIDE OUR OBSERVED FIELD, AND FOR NGC 5371, WE HAVE MEASURED THE r(CR) FOR THE SPIRAL ARMS, HENCE NO 

RATIOS CALCULATED. 



Object 


r(CR) 

II 


r(CR) 
kpc 


r(bar)e 
kpc 


r(CR)/r(bar) e 


kpc 


r(CR,)/r(bar) F 


IC 342 


344+2^ 

J —79 


6 5+?-S 

u - u -1.5 


5.8 


1 l+n'5, 
-0.3 


7.2 


9+°'i 


NGC 2403 


390 + L 


6 2 + I » 

u ' — 0.5 


7.8 


8 + I , 

w * — 0.1 


6.6 


9 + n , 
— o.l 


NGC 4294 


29+i 


1 9 + ^ 
± - a — 0.4 


1.7 


1 l + no 
± ' ± — 0.2 


2.5 


u -° — 0.2 


NGC 4519 






4.7 




4.7 




NGC 5371 






13.8 




13.0 




NGC 5921 




S.ol?;* 


6.3 


1 o+O.l 

J - J -o.2 


11.0 


7+ ' 1 


NGC 5964 


oo_ 13 




7.6 


u.y_ 2 


4.3 


1 6+ 10 
J -°-0.3 


NGC 6946 


254+4 


8.0±l Q 


7.6 


1 l +? 
J - l -o.3 


7.0 


1 l +? 


NGC 7479 




i3.o±i : g 


11.8 


1 1+ 01 


10.0 


i a+0.1 
0.4 


NGC 7741 


109±i 7 


5 5+ ? 


5.6 


l-Oto.2 


7.5 


7 +? 
u -'-0.2 



TABLE 6 

THE PATTERN SPEED C p , r(CR), AND ILR RADIUS FOR OUR SAMPLE GALAXIES. 



Object 




r{CR) 


Resonance radius 


£~2p 


r(OR) 


Resonance radius 


Structure 




km s — 1 arcsec - 1 


II 


// 


km s™ 1 kpc - 


kpc 


kpc 




IC 342 


n c;q+0.09 
U - OJ -0.09 


344±?| 




3ll 5 5 


6.5i;-j 




Bar 


NGC 2403 


34+° 14 
u -° -0.02 


390I32 




22+.? 






Bar 


NGC 4294 


2 92+ ' 14 
-0.81 


29l 4 6 




43+3 
* J -12 






Bar 


NGC 4519 


1 46+ - 44 


> 68 


20 (ILR) 


20l 6 8 


> 5 


1.5 (ILR) 


Spiral 


NGC 5371 


9 t-y+0.92 
z -°'-0.18 


93ii 6 




14±f 


17.0t ? 3. o 




Spiral 


NGC 5921 


1 59+ 024 
l.o»_ 24 






13+J 


8.ot°;§ 




Spiral/Bar 


NGC 5964 


o 0Q +0 .12 


trc:+38 
°°-13 




251^ 


7 o+ 4 - 5 

'•"-1.6 




Spiral/Bar 


NGC 6946 


n yq+o.19 

u -' J -0.19 




54 (OILR) 


25l 6 6 


8.0+i.o 


1.7(OILR) 


Spiral/Oval 


NGC 7479 


9 OO+0.47 
z.oo_ 47 


94i 6 25 




18l 3 3 


13-0±i;g 




Bar 


NGC 7741 


i ,,+0.48 
1 - lo -0.36 


109i' 7 


22 (ILR) 


19j* 


5 5+ ? 

°- J -i.o 


1.3 (ILR) 


Bar 



In the rightmost column, we specify the morphological 

is a measure of the intensity of the bar divided by 
the intensity of the inter-bar. As these authors have 
demonstrated, the contrast parameter is smaller in 
the region of galaxies where the bulge component 
dominates, then increases and reaches its peak at 
the middle of the bar, and falls at the end of the 
bar. 

3. The phase shift of the second Fourier c omponent 
(m 2 ) changes sign ( Zhang fc Buta||2007 1 . 



We note in Fig. [8] that not all these three criteria are 
fulfilled in our profiles, most likely due to presence of 
strong star formation knots, foreground stars, and dust 
lanes. However, we note that in some cases, in the 
proximity of r(bar) e , the first two criteria could hold, 



structure to which the pattern speed has been attributed. 

and in all cases, the phase of the second Fourier com- 
ponent crosses the zero line. We thus define the phase 
shift crossing radius as the r(bar)p- Table [5] shows that 
this method yields results sim ilar to the eflipticity pro- 
file analysis outlined in section [5. 1| and overall we derive 
comparable r(CR)/r(bar) ratios using either of the pho- 
tometric analysis methods. One prominent exception is 
the galaxy, NGC 4294, where the Fourier decomposition 
method unveils an extended 4-6 kpc oval structure which 
is not seen in ellipticity profile analysis or directly in the 
images. This oval structure could cause the severe under- 
estimation of our r(C R) / r{bar) p for this galaxy. 



5.3. Comparison With Simulations 
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From left to right: Projected ellipticity and position angle </> profiles for our sample galaxies, deprojected 2 nd (black) and 



4 (blue) Fourier component profiles, contrast parameter, and phase shift of the second component. In each panel, the vertical solid line 
indicates r(C-R), the diamond, bar radius derived from the ellipticity profile r(bar) e , and the vertical dashed line, the bar radius from the 
Fourier decomposition r(bar)p. In the <f> column, the horizontal dashed line shows the outer disk <f> from the near infrared image, and the 
dotted (red) line shows the p used when applying the Tremaine- Weinberg method. 
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Numerical simulations of barred galaxies predict the 
slowdown of bars due to the exchange of angular mo- 
mentum be tween the dark matter halo and the galactic 
disk ( e.g., Debattista fc Sellwood||1998| Athanassoula 
20031 IMartinez-Valpuesta et al.l|2006| |Sellwood fc De- 
battista||2006[ |Weinberg fe Katz|[2007| . Bars in the fast 
rotator regime are expected to have r(CR) /r(bar) be- 
tween 1.0 and 1 .4, except for simulations of [V alcnzucla 
fcKlypin| ( [2003l ) and |Tiret fc Combes| ( [2007l|2008p , which 
obtained slow rotators with r(CR) /r(bar) > 1.4. This 
ratio strongly relies on the measurement of the bar ra- 
dius, which is far from trivial to derive homogeneously 
from observations, however, the consistency check pre- 
sented in section [5] demonstrates that our morpholog- 
ically derived r(CR)/r(bar) ratios are consistent with 
those predicted by numerical simulations. The average 
value for our measurements is r(CR)/r(bar) = 1.1. As 
our main goal is to derive the f2 p , the r(CR)/r(bar) can 
serve as a useful inspection since its implication is that 
the bars in our sample late type spir al galaxies are fast 



bars, as predicted by models (e.g., Laine et al. 1998 
Debattista fc Sellwood||2000| |Rautiainen et al.||2008p . 

6. KINEMATIC RESULTS AND DISCUSSION 

6.1. Velocity Dispersion in Disks of Late-type Spiral 
Galaxies 

The galactocentric a profiles, illustrated in Fig.[l] show 
that for most galaxies a is almost constant through- 
out the observed field with some general hint of de- 
cline toward the 25th magnitude radius. NGC 5371 and 
NGC 7479 both show smaller Ha a values in the central 
ps 10" with the outer parts following the same trend as 
the rest of the sample. NGC 5921 has a deficit of ionised 
gas in the bar region, outside which it displays a very 
flat cr-profile. Consistent with the statistical results of 



Erwin & Sparke ([2002) for early-type bars and assum- 
ing that epicyclic approximation can be used here, we 
find that three of our 10 late-type spirals (NGC 4519, 
NGC 6946, NGC 7741) show evid ence for an inner Lind- 
blad Resonance ILR (see section 6.2 and Fig. [lj. In all 
three cases, the circular component of the observed V\ os 
(see section 3) shows a clear presence of a more rapidly 
rotating component within the ILR. Inside the ILR, the 
a values are generally higher than outside this region, 
with a particularly steep rise toward the nucleus, and as 
displayed in the rotation curves, the higher a values ap- 
pear to accompany more rapid rotation. This might be 
linked to the scenario where pseudo-bulges are caused by 
bar dissolution ( see the detailed review by Kormendy & 
Kennicutt|[2004 i. 

The results for our sample compare well with two- 
dimen s ional kinematic map s from Falco n-Barroso et al. 



(|2006[ ), rPeletier et al.| ( |2007[ ), and |Boker et al 
first order, our spatial resolution allows us to distinguish 
between local effects such as winds from star-forming re- 
gions and supernovae, and global secular evolution ef- 
fects. We find that the a values are h igher at the location 
of the star-forming Hn regions (c.f., |Allard et al^j [2006), 
as compared with regions where no strong star formation 
occurs. This result is similar to that which we found in 
Paper I, i.e. that locally outflows from Hn regions lead 
to a higher measured velocity dispersion close to the re- 
gions, while globally there is very little radial dependence 
of sigma, implying large scale uniformity throughout the 



disk in the summed inputs of the energy sources to the 
gas. 

6.2. Notes on Individual Galaxies 

IC 342 has angular size of more than 20 arcminutes and 
has a weak primary bar with r(CR) ss 4.7', in approxi- 
mately north-south direction ( Crosthwaite|2002" I . IC 342 
contains extremely strong molecular gas in the innermost 
parts of the bar, where the presence of an ILR (w 2") has 



been suggested (|Sheth et al.||2002| |Turner fc Hurt||1992 
Boker et al. 1997] T At this radius a star-forming rmg- 
like structure has been observed to include a « 1 5 Myr 
old superg iant-dominated central stellar cluster (Boker 
et al.|1997| . The near infrared spectra show no evidence 
for multiple stellar content within the bar, in support of 
this system being young. Although the exact molecular 
content in th e bar has not been est ablis hed, the CO ob- 
servat ions by Turner & Hurt ( 1992 ) and |Sakamoto et al 
dl999| suggest gas inflow due to the bar feeding the cir- 
cumnuclear starburst. Inside the ILR ring, a deficit of 
molecular gas content has been observed, and our obser- 
vations show a very weak Ha signal. We find an over- 
all patchy ionised gas distribution with a steeply rising 
rotation curve, fully consistent with both CO and Hi 
rotation curves which stay around 200 km s -1 out to 15' 



jRogstad et al.|1973l |Young 
|Crosthwaite et al.||2000| ). 



fc Scoville|1982l|Sofue|1996 



The bar f} D h a s been estimated to 40 km s 1 kpc 1 by 
Turner & Hurt|(1992l. Radial streaming motions above 
5 km s along one arm were detected by 



Newton 



fli980 l 

who also found r(CR) between 16 and 21 kpc, which 
they found consistent with a spiral density wave with Cl p 
kpc" 



= 10 km s 
the bar 



21 km s 



Crosthwaite et al. 
r 



(2000) derived 



p — *x iY±i± o kpc - based on morphology 
of the spiral structure and r(CR) at around 400" with 
no sign of an ILR. Applying the Tremaine- Weinberg 
method on the Ha emitting ionised gas, we derive 
n p = 31±1 km s" 1 kpc -1 and r(CR) « 6.5 kpc. We 
cannot verify the presence of the location of the ILR 
due to the large spatial bins which we use in order to 
increase the signal in the central regions (see Fig. [TJ. 
Finally, we confirm that the GHaFAS observations 
deliver a (V) / (X) value comparable to that when using 
the inner « 2' section of the FaNTOmM data. 

NGC 2403 is a high surface brightness member of the 
M 81 group with flocculent arms and bright H II regions 
and a tota l star formation rate « 0.1 M yr _1 (Kenni- 
cutt| [T989). In the inner kpc, this galaxy contains about 



half as muc h molecular gas as its ato mic hydrogen con 



tent (e.g., Thornley & Wilson 1995), suggesting that 



20081. To the formation of molecular hydroge n is constrained b y 



the low surface density of the gas (Heifer et al. 20031 



The central gas surface de nsities are below t h e threshold 



star formation density of |Kennicutt] ( |1989[ ) . Fraternali 
et al. (20011 have detected anomalous neutral hydrogen 



extending beyond the disk plane of the galaxy interpreted 
as galactic fountains powered by stellar winds from mas- 
sive stars and supernova explosions, which in turn could 
be triggered by the bar. Confirmed by a large sample of 
H II regions, there is a gradien t in O/H across NGC 2403 
with a slope of 0.1 de x kpc -1 ([Garnett et al .|l997[ ), sup- 
porting the study by Beckman et al. "( 1987p wEo found 
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that non-circular motions are not very efficient in mix- 
i ng th e ISM in this galaxy. Although |Schoenmakers et al.| 
(1 19971) found, using H I kinematics, that NGC 2403 is ax- 



isymmetric to a high degree, we find that the ellipse fits 
and Fourier analysis of the near infrared images reveal 
the presence of a sw 6 kpc bar ( see Fig. 

Our observed Ha rotation curve shows a moderate 
increase to w 60 km s" 1 in the central kpc followed 
by almost a solid-body type increase out to around 10 
kpc radius. The curve is con sistent with the Hi curve 
presented in Shostak (19731; however, due to strong 
feed-back from the numerous H II regions in the central 
w 200 pc radius (confirmed by the large a values), we 
cannot constra in the rotation curve within this region. 
Shostak (1973) derived the spiral pattern speed 10.5 
km s _ kpc - , whereas, using the Tremaine- Weinberg 
method, we derive the bar pattern speed f2 p = 22+^ 

km s" 1 kpc" 1 , and r(CR) at 6 kpc, with no sign of 
Lindblad resonances (Fig.nl). 

NGC 4294 appears to interact with NGC 4299 with a 
projected separation of 27 kpc. Both galaxies are sim- 
ilar in morphology, mass, luminosity, gas content, and 
both exhibit a high global star formati on rate (e.g., 
Koopmann fc Kenney 20041. Moreover, Chung et al. 
Q2007p found that due to its position in the Virgo cluster 
(0.7 JVIpc from M 87) interaction with the intra-cluster 
medium could be plausible. NGC 4294 hosts many Hn 
regions, a quiescent nucleus, and two rather faint floccu- 
lent spiral arms. The outer regions of the spiral struc- 
ture are consistent with enhanced star formation at the 
radius at which the gas surface densit y in the disk be- 
comes super-critical for st ar formation (Thornley & Wil- 



son 



l99^|Thornleyl[l996 l 



The observed Ha vel ocities show streamin g motions 

and this 

work) . 



along the spiral arms (Chemin et al. 2006 

Although the receding part of the disk appears 



slightly perturbed, the overall appearance of the velocity 
field is quite regular. Th e Ha rotation curve is consistent 
with that presented by Rubin et al. (19991. Using the 



rip — 



43 



+ 3 



Tremaine- Weinberg method, we derive 

km s" 1 kpc" 1 , with r(CR) 1.9 kpc in agreement 
with the location where the spiral arms emerge from 
the ends of the bar. Within the r(CR), the velocity 
dispersion follows a plateau, whereas outside this radius 
it decreases steadily. 

NGC 4519 hosts many Hll regions in the bar as well 
as the two partly broken spiral arms. The central re- 



gions exhi bit significant molecular gas content (Boker 
et al. 2002| ), although the star formatio n efficiency is con- 
stant in "the inner two effective radii (Rownd & Young 
19991. Our velocity field displays a S" -shaped zero- 
velocity curve over the central 4 kpc radius region of 
the galaxy. The rotation curve rises steeply in the cen- 
tral few arcseconds, which indicates a rapidly rotating 
component within ss 30", and it continues to increase to 
130 km s" 1 at 6 kpc radius . This is consistent with the 
optical rotation curve from Rubin et al. J 1999|) an d the 
H I position- velocity diagram of Helou et al. 1 1984 1 
We derive fi 



20^ 



km s 1 kpc , which we assign 



to outer spiral arms. We use this fi p to locate an outer 
ILR radius at about 1.5 kpc (20") and an inner ILR ra- 



dius at around 200 pc. The outer ILR radius is, to within 
the erro rs, consistent with the H a disk scale length de- 
rived by [Koopm ann et ap fl2006| ) . 

Similar to |Fathi e t al. (2007b), we apply the Tremaine- 
Weinberg method on the pixels interior to the ILR, 
and find that the inner bar is decoupled from the outer 
oval, as its pattern speed is f2 p = 45 km s" 1 kpc" 1 . 
This Op yields that the corotation of the inner bar is 
located at ILR of the outer oval. The method we have 
applied here enables us to determine the secondary Q p 
only with an uncertainty of 50%, since to zero-th order, 
the Tremaine- Weinberg method applied to the central 
region has to assume that the outer bar or spiral arms 
are stationary. This is not a realistic scenario, so we use 
this higher pattern speed for the inner bar only as an 
indication that we can detect the phenomenon, and that 
the pattern speed is higher than that of the outer bar. 

NGC 5371 is a grand design, possib ly post starburst 



spiral galaxy with a LINER nucleus (Rush et al. 1993 



|Koornneef||1993| |Elfhag et al.||1996[ ). It hosts a promi 
nent 2 kpc stellar bar sa 45 degrees from the orientation 
of the outer disk major-axis with no sig nificant change 
when comparing B and ii-band images (Eskridge et al. 
|2002| , and displays a star form ation rate of 0.9 M Q yr _1 
( |Gonzalez Delgado et al. 19 97]) . The rotation curve was 
measured by Zasov fe Sil'chenko] ( 1987 1 from which they 
derived an exponential main disk scale length of w 8.7 
kpc, and a massive bulge {M^uige > Mdisk)- Although 
the outer parts of their rotation curve are not well deter- 
mined, the inner regions agree well with our Ha rotation 
curve. 

Our observations display a deficiency of Ha emission 
across the bar, in agreemen t with the Ha images from 
Gonzalez Delgado k, Perez (19971. The velocity field 



displays clear disk-like rotation, with various "wiggles 
in the zero- velocity curve, often interpreted as stream - 
ing motions along spiral arms (Emsellem et al. 20061. 
The rotation curve rises steadily in the central 30 , 
and reaches 250 km s" 1 at 15 kpc radius. We derive 
n p = 14±[ km s" 1 kpc" 1 and the r(CR) w 17 kpc, i.e., 
93". The deficiency of the Ha emission from the stellar 
bar, combined with its non-optimal orientation compared 
with the outer disk major-axis, indicate that we are prob- 
ing the J7p of the spiral arms, and not the well-known 2 
kpc bar. 

We note that by placing corotation at 17 kpc, we bring 
the inner 4:1 resonance close to the end of the symmetric 
part of the spiral struct ure. At this region , in agreement 
with the predictions by Patsis et al. ( 1997 1, a bifurcation 
of the northern arm is seen in Fig. [TJ Furthermore, 
the association of the end of the symmetric part of the 
spirals with the inner 4:1 resonance is supported by the 
radial variation of the ratio of t he amplitudes of the 
m = 4 to the m = 2 components (Grosbol 1985 Patsis 



1991) 



NGC 59 21 has a ve ry poor molecular and neutral gas 
content (Verter 19851 and relatively poor ionised gas in 
a w 8 — 12 kpc prominent stellar bar which is enclosed 
by sharp inner ring-like structure. The bar contains two 
straight dust lanes, and the spiral arms host many H H 
regions (Gonzalez Delgado & Perez 1997 Aguerri et al. 
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1998 Rautiainen et al. 20051, and form the ring at the 
bar ends, just before they become very open, broad, and 
diffuse. At the north end o f the bar, the arm b e comes 
brighter as it leaves the ring ( Eskridge et al.|2002 ). Mar- 



tin k Friedli (19971 found that the stellar bar has an axis 



ratio of 0.5 and is oriented 18 degrees in the north-east 
direction. These authors also derived an asymmetric star 
formation efficiency < 0.2 M Q yr^ 1 in the bar, and seven 



times higher value in the circumnuclear region ( Gonzalez 
|Delgado et : aTJl997| >. 

Like in JNGC5371, the Ha emission is so weak in the 
bar of this galaxy that we have to make use of large 
spatial bins to measure any kinematic information. 
The rotation curve cannot be constrained within the 
bar. The Tremaine- Weinberg method thus delivers 
the pattern speed of the prominent spiral structure 
(Q p — 13^2 km s" 1 kpc" 1 ), with the spiral pattern 
r(CR) w 8 kpc (or 65", in good agree ment with the 



surface density (e.g., that in NGC5055 is a factor of 6 
lower than the gas density in M51). To summarise, our 
results for NGC 5964 imply an underlying density wave 
with tt p = 25±5 km s" 1 kpc" 1 . 

NGC 6946 is a grand-design barred spiral galaxy with 
three main gravitational distortions: a large oval with 
radius R sa 4.5', a R « 1' primary stellar bar with 
ellipticity 0.15, and a R w 8" nuclear bar with ellipticity 
0.4, which is almost perpendicular to the main bar. 
We have presented a detaile d analysis of t he reso nance 



structure in this galaxy in Fathi et al. (2007b I, and 



confirm that the maps presented in this paper are fully 
consistent with our previous work. 



NGC 7479 is an extensively studied starburst 



Hawarden et al. 1986 



yuti) w 8 Kpc jor do" , m good agree ment witn tne [2001 1 with a ma ssive, 2.2 x 10 ivi© s uenai uai wim axii 
{CR) = 71" of |Rautiainen et al.|[2008| ), which i s just ratio of 0.25-0.4 ( |Burbidge et al.|1960] |Martin|1995||. Al 



2000 Lainc 



Aguerri et al 
10 Mm stellar bar wit' 



alaxy 



axis 



k Zhang] 2009 1 . This suggests that the bar sh ares the 



outside the 58" radius of the strong stellar bar (|Buta| though there is no evidence of recent interaction, |Laine| 



same pattern speed as the spiral arms (e.g., Tagger 



eTaLl[T987 l 



k Heller (1999) have found many kinematic and morpho- 
logical features consistent with a minor merger event for 
this galaxy, proposing the hypothesis that this galaxy 
is on its way to evolve toward an earlier Hubble type. 



NGC 5964 hosts a very prominent and elongated bar 
surrounded by diffuse, extende d, and fragmented spiral 



arms with no regular pattern (Elmegreen k Elmegreen 



19871. This flocculent galaxy is rich in neutral Hydrogen 

content (M m /M T ot = 82%) and hosts a patchy circum- 119981). The mean star formation efficiency in the bar, al- 



106" (Elmegreen & Elmegreen 1985a; IQuillen et al. 11995 


Beckman k Cepa 


1990| and il v to be between 2 


7 and 


100 km s" 1 kpc" 1 


(|Laine et al. 


1998 


del Rio k Cepa 



nuclear structure ( Giovanelli k Haynes|1983| |Boker et al 
2002 1 . The Ha rotation curve displays almost a solid 



though asymmetric, is 0.4 M© yr 1 an d increases to 0.63 



body rotation in the inner 10 kpc. Our analysis shows 
that fl p = 25^5 km s" 1 kpc" 1 is well outlined in Fig. [Tj 
with the location of the r(CR) at 6.5 kpc. 

NGC 5964 displays the largest misalignment (36 de- 
grees) between the <p p and (f>k, which could be caused 
by the projection effect since this galaxy is near face-on. 
Such misplacement can also be seen as an indication of 
strong streaming along the spiral arms and prominent 
bar. In Fig. [I] we show that the r(CR) agrees with the 
radius where the patchy spiral structure starts, indicat- 
ing t hat the bar and spi ral arms rotate at the same rate 
(see |Tagger et ai]|1987[ ). 
Theoretic a l wor k" such as that by Elmegre en fc| 



Mq y r in the circumnuclear regions (Martin k Friedli 



19971, indicating inflow along the bar. 
Our Ha rotation curve is consistent with that de- 



rived by Garrido et al. ( 2005 1 with a prominent and 
rapidly rotating component in the central 30". This 
curve, as well as the angular frequency curve, cannot 
be determined in the central kpc region due to the 
weak Ha emission, and hence we cannot confirm the 

presence of the ILR reported by Laine et al. (|1998|) 

-l 



find fl r 



18±1 km s" 



kpc" 1 , and r(CR) 



We 
13 kpc. 



Our r(CR ) is alm ost twice the 55" found by Puerari 
k Dottori ( 1997j, the r(CR) sug gested by numerical 
simulations of |Laine et al.| (11998ft , and the 58" found 



Thomasson ( 1993 1 have shown that flocculent spirals 
can form via swing amplification of star formation 
patches, and thus should have no pattern speed. 
Although their proposed mechanism explains the patchi- 
ness of these galaxies, it does not exclude the presence of 
density waves in the galaxies where swing amplification 
builds patchy spiral structure. Su pport for such a 
scena rio comes from simulations by ( Wada k Norman 
2001 1 a s well as observational studi es~5y " Thornley k 



fro m potential-den sity phase-shift analysis carried out 



Mundyl (|1997[ ); |Grosbol k Patsisj dl998l); |ElmegreerT 
et aU pOOSp jotably, in 1NGC 5055, |Thornley k Mundy 



(19971 could outline the spiral arms by a spline fit 
to the _KT-band image and show that the non-circular 
motions in their residual H I velocity field follow these 
arms. These authors combined the photometry with 
kinematics to derive the pattern speed for the underlying 
density wave at 35-40 km s" 1 kpc" 1 , and conclude that 
in some flocculent galaxies the underlying density wave 
is difficult to observe due to their overall lower gas 



by |Buta k Zhahg| p009[ ). Using the r{CR) found by 
these authors, our angular frequency curve implies 
fl p sa 23 km s" 1 kpc" 1 , i.e., slightly above the value 
presented here. This would also make the bar with 
r(CR)/r(bar)p ~ 0.9, to within the errors, comparable 
with the value we present in Table [5] 

NGC 7741 has a prominent bar with axis ratio 0.36, two 
short and flocculent spiral arms including a significant 
amount of diffuse ionised gas and man y Hli regions (e.g 



Duval et al. 1991 Block et al. 2004 Laurikaincn et al. 



2004[ ). This bulge-less gala xy has a mean star form ation 
efficiency of 0.1 M Q yr -1 dMartin k Friedli||l997|), with 



lency ot U.l JYLq yr~" (Martin 
very weak C O content in the bar ( |Braine et al.| |1993). 
Duval et al. (19911 found strong ionised gas flows in the 
bar region and a ratio of 30% between the bar mass and 
that of the inner disk inside the bar radius. 
Our rotation curve agre es with the low-resolution 



curve of Duval k Monnet ( 1985 1 and Garrido et al 
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(2002). The Tremaine- Weinberg method results in 

main bar with 



r(CR) 



19; 



km s kpc for 



the 

l 



6. 5 kpc , c.f., fi p =31 km s kpc from Duval 



& Monnet (19851. We find the location of an outer ILK 
■ 1.3 kpc, mside which the rotation curve suggests 
presence of a dynamically cold rapidly rotating 
We find also an inner ILR at ss 300 pc, 



at ; 
the 

component 

and no hint of a higher pattern speed inside the outer 
or the inner ILR. In NGC7741, like in NGC7479, the 
ILR radii given here are based on the assumption that 
epicyclic approximation can be applied, and it should 
be noted that in strong bars, such as in NGC7741, 
where epicyclic approximation breaks down, t he actual 
resonances m ay be in different locations (e.g 
Teuben||2004| ). 



Regan & 



7. CONCLUSIONS 

Here, we have shown that an intensity- velocity map in 
Ha of the type obtained using a Fabry-Perot interferom- 
eter, can be successfully used to derive the gas kinemat- 
ics of disk galaxies with coverage and precision sufficient 
to derive the f2 p for those where star formation is well 
spre ad across the galaxy. Although this technique is not 



new ([H ernandez et al. 2005b ; E msellem et al.|2006| |Fathi| 
et al.||2007b = | 7 its precision and the number of galaxies 
analysed here give our results considerable value. It is 
clear that the Ha emitting gas does not act as a linear 
probe for the galaxy's surface density. However, only a 
fraction of the pixels in our maps are attributed to the 
compact H II regions. The Tremaine- Weinberg method 
applied in this way gives us data points with high spa- 
tial and velocity resolution over a major fraction of the 
galaxy disk, and our results demonstrate that, with cau- 
tion, the Tremaine- Weinberg method can be applied to 
deliver reliable pattern speeds in disk galaxies. 

Although very si mple , the radial tests that we have 
described in section |4.1| have proven useful to get a first 
order estimation of a possible secondary pattern speed. 
Fig. [7] shows that the pattern speeds decline steadily as 
larger sections of the field are covered, and exactly how 
much of the disk needs to be covered to get reliable val- 
ues depends very strongly on the size of the bars and/or 
spiral arms. 

Assuming that the epicyclic approximation can be 
used, the pattern speeds can be used to identify reso- 
nances, such as the ILR, in our sample galaxies. The 
presence of an ILR does not alone imply a decoupling of 
the region inside it, however, a n ILR is necessary for a 



decou pling of nested bars (e.g., Englmaier & Shlosman 
20041. We do trust the presence of secondary pattern 



speeds based on a combination of the radial behaviour 
of the (V) I ' (x) I sin(i), the presence of the ILR, and the 
presence of morphologically distinct features such as an 
inner bar. In galaxies without an ILR, we do not claim 
the presence of a secondary pattern speed, thus, once cer- 
tain of the mere presence of a secondary pattern speed, 
its presence can be confirmed using the full slit coverage 
but only for slits falling in the region corresponding to 
the extent of the inner bar. The numeric value of the 
secondary pattern speed, on the other side, is far more 
uncertain. 

We confirm our previous results from |Fathi et al.| 



ondary pattern speeds are as large as 50%, we do not 
include the values here, however, we note that the de- 
tected secondary patterns are of the order of two to 
three times higher than the main pattern speeds. Similar 



by |Corsini et al. 1 


2003) 


sellem et al. (|2006 


I, and 



the Ha maps of 10 late- type spiral galaxies can be listed 
as follows: 

• By sectioning the observed velocity fields in galac- 
tocentric rings, we have been able to illustrate that 
(V)/(X) is always highest at the centre, with a 
radially asymptotic decrease as pixels from larger 
sections of the disks are used, and converging to- 
wards the fl p when the full disk is covered. 

• Using the red continuum as our weighting function, 
we have found that the derived pattern speeds do 
not differ from those obtained when the Ha surface 
brightness maps were used. 

• The Tremaine- Weinberg method applied along the 
photometric major-axis gives least scatter fits as 
compared with extracting the parameters along the 
kinematic major-axis. This further supports the 
presence of strong non-circular motions in barred 
spiral galaxies which could artificially induce mis- 
alignment between the photometric and kinematic 
position angles. Such non-circular motions are ca- 
pable of transporting gas from the outer parts of 
the disks towards the circumnuclear regions, where 
secondary bars or inner disks can form. 

• Comparing the kinematically derived corotation 
radii for the bars with the bar radii independently 
derived from the morphology, we find that our 
r{CR)/r(bar) are consistent with predictions from 
numerical models using static or live dark matter 
halos or modified gravity. 

In six galaxies (IC342, NGC2403, NGC4294, 
NGC4519, NGC6946, NGC7479, and NGC7741), we 
have directly derived the f2 p of the main bar, and in 
four galaxies (NGC4519, NGC5371, NGC5921, and 
NGC 5964) we have been able to derive the f2 p of the spi- 
ral arms. Dynamically, bars and spiral arms can in teract 
directly wh en they are corotating (e.g., NGC 1365; Lind 



blad 
whic 

arms ove rlap (e.g., |Tagger et al.||1987| |Masset fc Tag- 



1999), or by means of non-linear mode coupling in 
case the r (CR ) of the bar and the ILR of the spiral 



ger|1997|). In two galaxies (NGC 5921 and NGC 5964) we 



(2007b I and since the uncertainties in deriving the sec- 



have found that the r(CR) of the spiral structure roughly 
coincides with the r(CR) of the bar, which suggests a rich 
variety of resonant interaction that have been predicted 
by theory, and our derived Sl p s complement the values 
for 10 late- type barred spiral galaxies in excellent agree- 
ment with f2 p s from numerical simulations using a live 
dark matter halo or modified gravity. 

We use epicyclic approximation, as commonly used 
in models, to identify the location of the resonances in 
our sample galaxies. Three of our ten late- type spirals 
(NGC 4519, NGC 6946, NGC 7741) show evidence for an 
ILR, two of which suggest the presence of a secondary 
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pattern speed. In all three galaxies, the rotation curve 
shows a clear presence of a more rapidly rotating compo- 
nent within the ILR, suggesting that they harbour sub- 
stantial amounts of the interstellar medium in the central 
regions. In all three cases, the Ha velocity dispersion is 
higher inside the ILR with a steep rise toward the nu- 
cleus. The combination of these two observational fea- 
tures supports the idea that the thickening of the gaseous 
component in the central region could build a bulge-like 
component in late- type spirals (see also Paper I). We 
note however that in strong bars such as in NGC 7741, 
epicyclic approximation could break down, and conse- 
quently, the location of resonances would not be valid. 

The present paper is the second in a series in which we 
will pursue detailed dynamical analysis of late-type spiral 
galaxies. Our current results will be complemented with 
a full harmonic analysis of the velocity fields followed by 
detailed numerical simulations, and more detailed esti- 
mates of the systematic effects due to non-continuity in 
the Ha emitting gas by applying the Tremaine- Weinberg 
method to high resolution simulations. Furthermore, a 
combination of these information with a robust estima- 
tion of the bar and spiral strengths could help us to better 
understand the evolution of structure in spiral galaxies. 
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